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ABSTRACT 


/ 

— >>  The  steady,  two-dimensional  flow  of  an  Incompressible  fluid 
over  a  backward-facing  step  is  computed  by  a  finite  difference 
procedure.  The  flow  depends  upon  the  Reynolds  number.  Re,  based 
on  the  inlet  flow  conditions,  and  on  the  dimensionless  step 
height,  s.  Spatially  resolved,  accurate  solutions  are  obtained 
for  a  range  of  s  for  Reynolds  numbers  as  high  as  350.  The  pri¬ 
mary  flow  feature  is  a  steady  separated  region  immediately  in  back 
of  the  step.  Additionally,  a  secondary  separated  vortex  can 
appear  on  the  top,  straight  wall  of  the  channel.  A  region  is 
delineated  in  the  (s,Re)  plane  where  laminar  separation  occurs  on 
the  top  surface.  It  is  concluded  that  there  is  no  flow  reversal 
on  the  top  surface  if  Re  is  less  than  about  200,  regardless  of  the 
step  height  s.  The  occurrence  of  this  secondary  vortex  is 
associated  with  laminar  separation  of  the  top  boundary  layer  pro¬ 
duced  by  an  adverse  pressure  gradient  along  the  top  of  the 


1.  Introduction 

The  flow  field  of  an  Incompressible  fluid  over  a  backward -facing  step 
exhibits  many  of  the  nonlinear  features  which  are  Intrinsic  to  viscous  fluid 
motion.  Because  of  this  and  the  frequent  occurrence  of  a  sudden  expansion  in 
engineering  applications,  there  have  been  several  experimental  and  theoretical 
studies  of  this  flow  situation.  The  recent  paper  by  Armaly,  Durst,  Pereira  A 
Schonung  (1983)  Includes  a  comprehensive  literature  review  for  the  problem. 
Figure  1  is  a  sketch  of  the  geometry  and  the  possible  separated  flow  struc¬ 
ture.  The  step  height  is  s*  (the  superscript  *  refers  to  a  dimensional 

it  it  it 

quantity)  while  the  gap  size  Is  1  ;  thus  the  expansion  ratio  is  l:(l+s  /I  ). 
Flow  separation  occurs  at  the  tip  of  the  step,  and  the  primary  eddy  reattaches 
at  a  distance  xj*  from  the  Inlet.  This  distance  Is  a  function  of  the 
expansion  ratio  and  the  Reynolds  number  Re.  At  sufficiently  large  Re, 
secondary  separation  on  the  top  boundary  takes  place  (see  Figure  1)  at  x2* 
and  the  flow  reattaches  at  X3*.  These  phenomena  were  observed  by 
Armaly  et  al .  for  an  expansion  ratio  of  1:1.94  and  confirmed  by  the 
computation  of  Mansour,  Kim  A  Moln  (1983),  who  compute  the  flow  corresponding 
to  an  expansion  ratio  1:  2. 

While  the  major  flow  reversal  originates  due  to  the  sharp  edge  of  the 
step.  It  Is  evident  that  the  top  separation  comes  about,  at  least  In  part, 
because  of  the  adverse  pressure  gradient  associated  with  the  expanding  flow 
created  by  the  dividing  streamline  (which  reattaches  at  xj*).  In  addition, 
the  Reynolds  number  of  the  flow  has  to  be  sufficiently  large  to  bring  about  a 
boundary  layer  structure  with  an  associated  decrease  In  the  skin  friction, 
which  together  with  the  Increasing  pressure  leads  to  top  separation  at  x2* 
(cf.  Brown  A  Stewartson  1969).  However,  it  should  be  noted  that  the  adverse 


pressure  gradient  will  get  weaker  with  increasing  Re  as  xj  increases  with 
Re.  Thus,  for  a  given  expansion  ratio,  it  may  be  theorized  that  top  separa¬ 
tion,  if  any,  will  be  confined  to  a  finite  range  (above  zero)  of  Re. 

At  sufficiently  high  Re,  Armaly  et__aK  (1983)  observe  and  compute  yet 
another  flow  separation  and  reattachment  on  the  lower  surface  at  x4*  and  xg* 
(Figure  1).  The  discrepancy  between  their  experimental  and  numerical  results 
however  was  large.  This  is  due  to  the  large  inaccuracy  and  false  diffusion, 
associated  with  high  Re  computations,  of  the  upwind  scheme  they  used.  The 
complicated  picture  of  the  flow  field  sketched  in  Figure  1,  coupled  with  the 
fact  that  the  physical  model  extends  to  infinity  in  the  flow  direction, 
clearly  necessitate  the  use  of  a  highly  nonuniform  computational  mesh  to  ade¬ 
quately  resolve  the  different  regions  of  sharp  velocity  gradients.  Such  a  mesh 
can  be  designed  after  these  regions  have  been  identified  in  exploratory  compu¬ 
tations.  This  is  the  approach  we  take  in  this  work.  However,  it  is  recog¬ 
nized  that  some  form  of  an  adaptive  grid  refinement  may  be  an  attractive  and 
computationally  more  effective  alternative  (cf.  Berger  &  Oliger  1984). 

The  objectives  of  the  present  work  were  to  investigate  the  features  of 
the  flow  for  a  wide  range  of  s  than  is  presently  available,  with  particular 
attention  paid  to  the  causes  and  details  of  the  secondary  separation  on  the 
top  surface. 


2.  Mathematical  model  and  solution  technique 


The  nond 1 men slonal  Navler-Stokes  equations  are 


v  »  0, 

*** 

(la) 

2 

(vv)  *  -7p  +  7  v/Re. 

(lb) 

The  length  and  velocity  scales  are,  respectively,  1*  and  U,  where  U  is 
the  average  flow  velocity  at  the  inlet.  Thus  the  Reynolds  number  Re  is  Ul*/v 
where  v  is  the  kinematic  viscosity.  Note  that  this  definition  of  Re  is 
different  from  that  of  Armaly  et  al..  In  the  usual  two-dimensional  Cartesian 


notation,  the  boundary  conditions  are 

!(0,0)  ;  0  <  y  <  s, 

(2a) 

(6(y-s)( 1+s-y) ,0) ;  s  <  y  <  1+s, 

3X  v(-,y)  *  0,  (2b) 

v(x.O)  =  0,  (2c) 

v(x,l+s)  -  0.  (2d) 


Here  s  is  the  nondimensional  step  length  (s*/l*)  and  we  have  assumed  a  fully 
developed  parabolic  velocity  profile  at  the  inlet. 

Finite  difference  solutions  of  equations  (1-2)  are  computed  on  the  stag¬ 
gered  grid  shown  in  Figure  2.  The  computational  region  is  divided  into 


rectangular  control  volumes  with  the  grid  points  located  at  the  geometric  cen¬ 
ters  of  these  cells.  Additional  boundary  points  are  Included  where  the 
boundary  conditions  (2)  are  imposed.  The  number  of  the  x  grid  points  (control 
volumes)  are  NX(NX-2)  of  which  MX  control  volumes  are  allocated  to  a 
distance  -  x^.  Similarly,  the  number  of  y  grid  points  (control  volumes) 
are  NY(NY-2)  of  which  MY  control  volumes  are  allocated  to  the  gap.  As  shown 
in  Figure  2,  there  are  two  uniform  grids  In  the  y  direction;  the  grid  covering 
the  primary  separated  region  is  also  uniform  to  a  distance  x^.  The  grid 
increases  gradually  In  the  x  direction  (for  a  total  of  NX-2-MX  control 
volumes)  to  a  length  nxL.  We  set  n  at  3  which  was  found  to  be  sufficient  by 
Armaly  et_  aj_.  (1983)  and  Mansour  et_£k  (1983). 

The  details  of  the  difference  scheme  and  the  iteration  procedure  we  use 
to  solve  the  difference  equations  are  described  by  Patankar  (1981).  The 
finite  difference  equations  are  obtained  by  Integrating  equations  (1)  over  the 
control  volumes  with  assumed  local  linear  variations  in  all  of  the  primitive 
variables.  The  convective  and  diffusive  fluxes  are  approximated  by  a  power- 
law  scheme.  A  llne-by-llne  iteration  to  solve  the  discretized  equations  is 
used  with  one  iteration  comprising  four  double  sweeps  of  the  field.  Under¬ 
relaxation  in  solving  for  u  and  v  was  required;  a  relaxation  factor  of 
0.85  was  used  throughout.  Solutions  were  constructed  on  a  number  of  non- 
uniform  grids.  In  Table  1  we  list  the  parameters  of  8  such  grids  together 
with  the  VAX  11/780  computer  time  requirement  for  one  iteration.  Convergence 
was  assumed  when  the  maximum  change  in  u  or  v  over  one  Iteration  was  less 
than  10"^.  About  300  iterations  were  needed  for  convergence. 


3.  Numerical  results 

a.  Numerical  accuracy 

Solutions  of  equs.  (1-2)  were  found  for  a  range  of  (s,Re).  Before  a 
specific  run,  an  estimate  of  x^  was  made  In  order  to  ensure  a  sufficiently 
long  computational  region.  In  Table  II  we  list  some  of  the  important  details 
of  the  flow  fields  and  the  particular  grid  on  which  that  solution  is 
obtained.  In  Table  II  4*  is  the  stream  function  and  x^,  X2  and  X3  are  the 
values  of  x  where  u,  computed  on  the  grid  line  just  adjacent  to  a  wall, 
changes  sign.  In  order  to  be  consistent  with  previous  work,  we  have  expressed 
these  distances  normalized  with  respect  to  the  step  height,  s. 

A  general  observation  from  Table  II  Is  that  finer  grids  than  those  used 
here  are  needed  to  produce  accurate  solutions  at  the  larger  values  of  s  or 
Re.  Indeed,  the  length  of  the  primary  bubble  x^  Increases  with  Re  at  any 
value  of  s,  at  sufficiently  low  values  of  Re,  as  expected.  Armaly  et  al . 
(1983)  report  that  for  s  ■  0.94,  Xj  begins  to  decrease  as  Re  Increases  above 
200.  Mansour  et  al.  (1983),  for  s  *  1,  show  that  this  is  due  to  insufficient 
numerical  resolution.  Figure  3  shows  the  dependence  of  Xj  on  Re.  We  have 
added  our  results  to  those  of  previous  workers,  who  present  results  only  in 
graphical  form.  It  Is  evident  that  we  need  a  finer  mesh  than  M7  to  resolve 
flow  fields  above  Re  *  250  for  this  range  of  s.  It  can  also  be  learned  from 
Table  II  that  if  a  decrease  of  Xj  with  Re  occurs  with  results  from  some  grid, 
results  with  finer  grids  show  that  the  reverse  Is  correct.  Thus  the  decrease 
of  xj  with  Re  reported  by  Armaly  et_£l_.  Is  clearly  due  to  Insufficient 
resolution.  It  is  also  Interesting  that  our  80x50  nonuniform  mesh  apparently 
gives  results  for  Re  <_  250  which  are  comparable  to  those  of  Mansour  et_^l_.  who 
used  a  130*130  uniform  mesh. 


Separation  on  the  top  surface,  should  It  occur,  may  be  missed  If  the  grid 
is  not  sufficiently  fine  in  the  y  direction.  This  is  clear  from  the  results 
listed  for  s  *  0.94  and  Re  *  200  where  results  with  M3  (30  y-grid  points)  show 
no  flow  reversal  while  those  with  M4  (40  y-grid  points)  indicate  top  boundary 
layer  separation.  Even  though  the  results  are  not  as  accurate  as  we  would 
like  to  report,  it  is  quite  clear  that  we  can  state  whether  the  top  bubble 
exists  but  cannot  report  more  than  2  place  accuracy  for  X£  and  X3. 

In  summary,  the  details  of  the  primary  vortex  are  resolved  to  an  accuracy 
of  1-2%,  and  those  of  the  secondary  vortex  to  approximately  10%. 

b.  Primary  Vortex 

Figure  3  gives  our  main  results  on  the  characteristics  of  the 
primary  vortex,  which  for  the  range  of  s  covered  in  this  study,  is  of  a 
similar  nature  to  that  reported  by  Armaly  et  al .  and  Mansour  et_  al_.  Our 
computations  do  not  extend  to  a  sufficiently  high  Reynolds  number  to  suggest 
the  nature  of  the  eddy  growth  as  Re  -*■  «.  While  the  experiments  of 
Armaly  et_^K  and  theoretical  models  (see  Acrivos  A  Schrader  (1982)  for  a 
discussion)  suggest  that  ~  Re  as  Re  ♦  »,  no  numerical  results  to  date  are 
sufficiently  accurate  at  high  Re  to  allow  the  question  to  be  unequivocally 
answered.  We  do  note,  however,  that  since  our  mesh  M4  utilizes  a  resolution 
in  the  primary  vortex  of  comparable  quality  to  the  130x130  uniform  mesh  of 
Monsour  et_^l_.,  and  since  M4  is  of  Insufficient  accuracy  for  Re  >  250,  the 
decrease  in  the  rate  of  growth  of  xj  observed  by  the  latter  authors  is  of 
questionable  validity. 


c.  Secondary  vortex 

Armaly  et_£l_.  were  evidently  the  first  workers  to  conduct  a  careful 
study  of  the  secondary  vortex.  However  their  results  are  for  a  single  value 
of  s,  and  we  became  interested  in  the  dependence  of  the  nature  and  occurrence 
of  secondary  separation  on  this  parameter.  Our  main  results  are  shown  in 
Figure  4,  where  we  plot  the  approximate  locus  of  the  curve  in  the  Re,s  plane 
separating  the  regions  in  which  the  secondary  vortex  does  and  does  not 
occur.  An  interesting  feature  of  this  diagram  is  the  double-nosed  nature  of 
the  curve  near  s  *  1.0. 

This  figure  may  be  understood  by  first  interpreting  the  limits  s  *  0,  •. 
s  *  0  corresponds  to  channel  flow  which  evidently  has  no  steady  separated 
regions,  regardless  of  the  Reynolds  number.  Thus  the  locus  cannot  intercept 
the  line  s  *  0.0,  although  it  possibly  can  become  tangent  to  it  as  Re  ♦  «. 

The  other  limit,  s  ♦  -,  corresponds  to  a  wall  jet  issuing  into  stagnant 
surroundings,  and  again,  we  have  no  evidence  that  steady  separated  flows  exist 
in  that  case  as  well. 

As  we  argued  in  the  introduction,  the  occurrence  of  the  secondary  vortex 
Is  most  probably  a  boundary  layer  phenomenon  which  will  occur,  for  finite  s, 
for  a  sufficiently  large  Reynolds  number.  The  pressure  field  associated  with 
the  deceleration  of  the  expanding  flow  between  the  top  wall  and  the  zero 
streamline  connecting  Xj  with  the  lip  of  the  step  will,  at  high  Re,  contain  a 
region  of  adverse  pressure  gradient.  As  the  decelerating  boundary  layer  along 
the  top  wall  is  acted  upon  by  viscosity,  it  will  eventually  separate,  causing 
a  flow  reversal  and  a  diminishing  of  the  adverse  pressure  gradient  due  to  the 
growth  of  the  secondary  vortex. 
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All  of  our  available  numerical  evidence  is  in  support  of  this  interpreta¬ 
tion.  First,  we  note  that  in  all  cases,  the  point  of  separation  xg  is 
upstream  of  the  reattachment  point,  xj,  of  the  primary  vortex.  This  is 
consistent  with  the  above  argument  since,  for  x  »  x^,  the  pressure  gradient 
becomes  favorable  as  the  flow  relaxes  to  a  plane-Poiseui 1 le  flow.  Second,  we 
can  examine  the  pressure  profile  on  the  top  wall.  Shown  in  Figure  5  is  a 
series  of  profiles  showing  the  pressure  gradient  along  the  top  wall  as  a  func¬ 
tion  of  downstream  distance.  These  particular  profiles  are  for  s  =  0.8,  but 
other  values  of  s  give  analogous  results.  As  can  be  seen,  the  pressure  gra¬ 
dient  becomes  unfavorable  (positive)  a  short  distance  from  the  step.  As  the 
Reynolds  number  increases,  the  length  of  the  region  of  unfavorable  gradient 
increases  until  laminar  separation  occurs  on  the  top  face.  This  is 
accompanied  by  the  occurrence  of  two  inflection  points  in  the  pressure 
gradient.  Thus  comparing  Figures  4  and  5,  we  see  that  for  s  =  0.8 
Re  *  275  is  the  first  case  for  which  two  inflection  points  occur.  Above 
Re  =  275,  the  region  of  decreasing  pressure  gradient  grows,  as  does  the  length 
of  the  region  between  inflection  points.  This  latter  behavior  is  due  to  the 
growth  of  the  secondary  vortex.  These  points  are  emphasized  in  Figure  5,  where 
for  Re  =  325  the  quantities  xj,  X2»  X3  are  marked  on  the  curve.  The  first 
inflection  occurs  at  x  *  ,  while  the  second  occurs  slightly  upstream  of 

X3.  The  strong  correlation  between  the  pressure  gradient  and  the  appearance 
and  disappearance  of  the  secondary  vortex  is  further  illustrated  by  examining 
the  pressure  gradient  for  Re  =  225  as  a  function  of  s.  As  can  be  seen  from 
Figure  4,  such  a  series  of  profiles  traverses  regions  of  separated  and 
unseparated  flow  twice  as  s  is  increased  from  0.85  to  1.10.  The  family  of 
pressure  gradient  curves  is  shown  in  Figure  6,  and  as  can  be  seen,  there  is  an 
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excellent  correlation  between  the  double  inflection  points  in  the  pressure 
gradient  and  the  occurrence  of  top  separation. 

From  these  studies,  we  may  safely  conclude  that  the  occurrence  of  the 
secondary  vortex  is  a  laminar  separation  of  the  boundary  layer  on  the  top  wall 
due  to  the  adverse  pressure  gradient  set  up  by  the  deceleration  of  the  forward 
moving  fluid,  which  in  turn  is  due  to  the  growth  of  the  primary  vortex  with 
Reynolds  number. 

Some  of  the  detailed  features  of  Figure  4  remain  unexplained  to  us.  The 
occurrence  of  the  nose  at  s  *  1.05  may  be  partially  understood  by  examination 
of  the  dependence  of  xj  on  s  at  a  fixed  Reynolds  number.  Figure  7  shows  a 
number  of  such  curves,  and  as  can  be  seen,  in  all  cases  in  which  there  is  to 
separation  (  Re  >  200  ),  x^  Is  a  local  minimum  for  s  *  1.05.  (The  same 
statement  holds  true  If  the  quantity  Xj/s.as  opposed  to  Xj,  Is  examined.) 

This  Indicates  that  the  unfavorable  pressure  would  be  expected  to  be  a  maximum 
for  s  *  1.05,  leading  to  separation  at  the  lowest  value  of  Re  for  which  it 
occurs.  The  other  features  of  Figure  4  such  as  the  nose  at  Xj  *  0.94,  would 
appear  to  depend  on  more  subtle  details  of  the  flow. 

Finally,  we  note  that  for  s  ~  0(1),  the  secondary  vortex  cannot  grow  to 
be  very  deep,  although  It  apparently  can  grow  to  be  quite  long.  The  reason 
for  this  is  simply  that  should  It  penetrate  any  substantial  depth  into  the 
channel,  it  would  modify  and  possibly  eliminate  the  adverse  pressure  gradient 
which  Is  responsible  for  Its  occurrence.  By  similar  reasoning,  the  tertiary 
vortex  observed  by  Armaly  et  al .  Is  probably  a  turbulent  phenomena,  at  least 
for  s  ~  0(1),  since  the  secondary  vortex  never  grows  deep  enough  to  cause  a 
secondary  region  of  decelerating  flow,  with  another  associated  adverse 
pressure  gradient. 
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Table  I 

Mesh  Definition 


Mesh 

NX 

NY 

MX 

MY 

Ml 

30 

20 

15 

9 

M2 

40 

30 

20 

14 

M3 

50 

30 

25 

14 

M4 

60 

40 

30 

19 

M5 

70 

40 

35 

19 

M6 

80 

40 

40 

19 

M6* 

80 

40 

35 

19 

M7 

80 

50 

40 

24 

in  VAX  11/780  (UNIX)  seconds  per  iteration. 


CPU1 


3.49 

7.31 

9.30 

15.34 

18.03 

21.41 

21.41 

32.80 


12 


s 

0.7 


0.75 


0.80 


Table  II 

Variation  of  results  with  computational  mesh 


Re 

XL 

Mesh 

Xj/S 

x2/s 

x3/s 

300 

7.5 

M4 

4.37 

10.11 

325 

7.5 

M4 

4.36 

10.61 

10.01 

11.58 

M5 

4.34 

10.66 

10.05 

11.80 

M6 

4.32 

10.73 

10.18 

11.86 

350 

8.5 

M4 

4.37 

10.92 

9.52 

13.08 

M5 

4.35 

11.04 

9.68 

13.25 

M6 

4.33 

11.12 

9.81 

13.35 

275 

7.5 

M4 

4.80 

9.69 

300 

7.5 

M4 

4.79 

10.21 

9.35 

11.51 

M5 

4.77 

10.28 

9.40 

11.69 

MS 

4.74 

10.32 

9.50 

11.77 

325 

8.0 

M4 

4.79 

10.55 

9.04 

12.95 

M5 

4.77 

10.88 

9.19 

13.10 

MS 

4.75 

10.74 
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FIGURE  CAPTIONS 


Figure  1.  Sketch  of  the  flow  field  with  regions  of  flow 
reversal  identified. 


Figure  2.  The  computational  nonuniform,  staggered  grid.  There 
are  NX,  NY  mesh  points  in  the  x,y  directions, 
respectively. 


Figure  3.  The  length  of  the  primary  vortex,  x^/s,  as  a 
function  of  Reynolds  number. 


Figure  4.  Onset  of  laminar  separation  on  the  top  surface.  The 
crosses  stand  for  pairs  of  (s,Re)  where  there  is  no 
flow  reversal  while  the  circles  denote  separation. 


Figure  5.  Pressure  gradient  as  a  function  of  downstream 
distance:  s  *  0.8  with  Re  as  a  parameter.  The 
curves  have  been  shifted  in  x  for  clarity. 


Figure  6.  Pressure  gradient  as  a  function  of  downstream 
distance:  Re  *  225  with  s  as  a  parameter. 


Figure  7.  Primary  separation  distance  Xi  as  a  function  of  step 
height  s,  with  Re  as  a  parameter.  The  error  bars 
denote  the  maximum  error  due  to  the  variation  of 
grid  sizes  among  the  data  plotted. 
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